Thermal conductivity across a twin boundary in d-wave superconductor 



in 

(N 



o 
o 

u 

Oh' 



X3 
O 

o 



> 

00 
(N 

a\ 
o 

On 
-I— > 



I 

O 

o 



X 
S3 



M. E. Zhitomirsky and M. B. Walker 
Department of Physics, University of Toronto, Toronto, Canada M5S 1A7 

(February 1, 2008) 



We consider excitation spectrum near a twin boundary in an orthorhombic d + s superconduc- 
tor. The low-energy spectrum is highly sensitive to the presence of the small amount of s-wave 
component. Robustness of the bound states at the Fermi level with respect to impurities and an 
extra boundary potential is investigated. The role of Andreev transmission process for the low- 
temperature thermal conductivity across twin boundaries is studied for an impure superconductors. 
At very low temperatures the bulk part and the twin boundaries part of the thermal conductivity 
have similar linear-T dependences, whereas at intermediate temperatures the two contributions 
behave like and T^, respectively. 



PACS numbers: 74.25. Fy, 74.72. Bk 



I. INTRODUCTION 

The bound and extended quasiparticle states at twin 
boundaries in an orthorhombically twinned d-wave su- 
perconductor (like, supposedly, YBa2Gu307_2;) have re- 
cently been studied for the clean case. El An areal density 
of low energy twin-boundary bound states proportional 
to the magnitude of the s-wave component of the super- 
conducting gap (which is a linear combination of tetrag- 
onal s-wave and d-wave componentsu'EI) was found, and a 
process of Andreev transmission (particle-hole conversion 
in transmission across the twin boundary) was shown to 
be the essential mechanism of transporting heat across 
the twin boundaries in the limit of sufficiently low tem- 
peratures. 

The present article considers the effect of disorder 
on the twin-boundary bound and extended states. A 
number of studies have shown the presence of impu- 
rities gives rise to a significant bulk density of zero 
energy excitations in the d-wave superconductor and 
that these hajiie|^^gniiicant effect on the bulk thermal 
conductivity.tau-ffltJ'LrEjO This paper shows that the en- 
hanced density of states at the Fermi level affects also 
the transport of heat across the twin boundaries. Fur- 
thermore, in addition to giving a finite lifetime to the 
extended bulk excitations, impurities transform the twin- 
boundary bound states into decaying resonances. 

Theoretical calculations (see, e.g., Ref. ^ show that 
for a d-wave superconductor the low temperature a-b 
plane penetration depth varies linearly with T in the 
clean limit, but that impurity scattering can change this 
variation to T^. Since this result depends essentially on 
the fact that the gap displays lines of nodes, and not on 
the d-wave symmetry, it is applicable to YBa2Cu307_a: 
(YBCO), where the orthorhombic Aig gap is considered 
to be a mixture of s-wave and d-wave components. In the 
Born limit, the amount of scattering required to produce 
the variations of the penetration depth is so large that 
Tc will be dramatically reduced, but in the unitary limit 



impurity scattering can give the temperature depen- 
dence without gr|eatly affecting Tc- Very clean supercon- 
ductors do showLJ the low temperature a-b plane pen- 
etratiorL-deoth varying linearly with T whereas detailed 
analysialj'Ej of the results from a number of other sam- 
ples reveals a T^-behavior. For these reasons, much re- 
cent work has analyzed the effect of the unitary-limit im- 
purity scattering on the properties of high-temperature 
superconductors. 

Since a somewhat intuitive approach will be used be- 
low for calculation the thermal conductivity of the twin 
boundaxies (e.g., in comparison with the Ambegaokar- 
Griffintj approach to the calculations of bulk thermal 
conductivity), we review here the relevant qualitative 
ideas and their application to the description of other 
properties of high-temperature superconductivity. For a 
clean d-wave superconductor the density of states for the 
excitations in the superconducting state varies linearly 
with excitation energy at low energies and can be writ- 
ten in the form 



NsiE)/No^E/Ao 



(1) 



where Nq is the normal state density of states and Aq 
is the maximum gap. In the strong scattering limit, the 
energy of low-lying excitations with wave-vector k ac- 
quires an imaginary part, i.e., becomes + iF, thus 
the excitation possessing a lifetime r = l/(2r) {h = 1). 
The quantity F is given by the solution of the following 
equation (see, e.g., Ref. ||) 



mAo 



2Aro ln(4Ao/F) ' 



(2) 



where is the impurity concentration. Note that this 
result for the damping of an excitation is independent of 
the wave-vector and energy, and is valid for quasiparti- 
cles with energies less than F. Since, by the uncertainty 
principle, the low-energy excitations will have a spread 
in energy of the order of F, excitations with i?k from zero 
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up to r should contribute to the density of states at zero 
energy. Thus, taking account of Eq. (|l|), we have 



II. QUASIPARTICLE STATES NEAR TWIN 
BOUNDARY 



NsiE)/No^T/Ao 



(3) 



for the density of superconducting excitations at zero en- 
ergy in the presexice of unitary-hmit impurity scattering. 
The exact resuhEl differs from this by a logarithmic fac- 
tor which varies weakly with impurity concentratioH and 
which is often neglected in qualitative discussions .H The 
low temperature specific heat associated with this con- 
stant density of states at low energies is (fc^ = 1) 



27r2 r 
Cv = -^^o-^T . 
6 Aq 



(4) 



These ideas allow an estimate of the bulk thermal con- 
ductivity to be made asM 



3 Ao 



T. (5) 



Note that, like the result for the electrical conductivity in 
the unitary limit of impurity scatteringjj this is a univer- 
sal result. The parameter F, which depends on n^, cancels 
out and the thermal conductivity is independent of the 
impurity concentration and the strength of impurity po- 
tential. Recently, the universal heat transport has been 
observedliS in untwinned crystals YBa2(Cui_a;Zn2:)306.9. 
Quantitative comparison of these data with the theory 
suggests that quasiparticle scattering on Zn impurities 
in YBCO is close, but does not coincide completely with 
the unitary limit. 

At higher energies F <C -E <C Aq the density of states 
reverts to it's clean- limit value (0), whereas the relax- 
ation rate of superconducting quasiparticles acquires en- 
ergy dependence Ts — TnAq/E (Ref. where the nor- 
mal state relaxation rate in the unitary limit is F^r = 
ni/nNo- In the corresponding temperature range F <C 
T ^ Ao the electronic thermal conductivity can be cak 
culated by an elementary Boltzman equation approachO 



y-v {EkVk cos gk)- 

k 



TVs 



9Ek 



(6) 



provided the correct energy dependence of Fs given 
above is taken into account. Using the fact that near 
the gap nodes the quasiparticle group velocity reduces to 
Di?, we find after integration 



7n^ v],No 

15 fa,a: 



(7) 



Thus, above the universal linear-T behavior at very low 
temperatures the heat current in d-wave superconductor 
varies as , which is consistent with the numerical data 
presented in Ref. M. 



In a tight-binding model on a square lattice (suppos- 
edly reproducing in a more or less correct fashion a single 
Cu02 plane in YBCO), the Green function describing the 
excitations in the superconducting state is determined by 
the equation 



eik] G{k,j, z) = I5i- 



(8) 



Singlet state pairing is assumed and a "hat" indicates a 
2x2 matrix in particle-hole space. Disorder is included in 
this model by considering impurities on randomly chosen 
lattice sites, each impurity being described by an addi- 
tional potential u. Impurity scattering is treated in the 
self-consistent t-matrix approximation. As a result, the 
quantity z must be determined from the equation 



where 



niu'^Qn{z) 
l^u^g^izf 



Qo{z) = \TiG{i,i,z) 



(9) 



(10) 



The right hand side of Eq. ( |lO| ) is assumed to be indepen- 
dent of site index i, z is the usual Matsubara frequency, 
and the s-wave component of the gap is assumed negligi- 
ble in comparison with the c?-wave component in arriving 
at Eq. (|). 

In the unitary limit {u oo) detailed study of Eq. (||) 
shows that for low energy excitations, i.e., excitations 
having an energy less than F, the quantity z is given by 



z±iT 



(11) 



where F has been defined in the Introduction, and the 
upper and the lower signs refer to the upper and lower 
half plane, respectively. 

The matrix Cik contains the same quantities that were 
used in our previous modeLof a clean superconductor 
containing a twin boundary:EI 



Cifc = T3 {-Uk - fJ-Sik + UoSikSiB) + fiAik 



(12) 



where SiB is unity for i on the twin boundary and zero 
otherwise and A^/j is assumed real. 

When the site indices are incorporated into the matrix 
notation, Eq. (||) becomes 



[iz - e] G = i . 



(13) 



where z is given by Eq. (|ll|). Since e is a real symmetric 
matrix, it can be diagonalized by a unitary transforma- 
tion, say S, giving 



G S\iz- A)-^S 



(14) 
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where A is the diagonal matrix of the eigenvalues of e. 
The excitation energies are found from the poles of G 
when z — > £' + iO. Thus, in the low-energy limit, if is 
an eigenvalue of e, there will be an excitation with energy 



E = A„; 



(15) 



Therefore, all low energy excitations, including those 
which would be bound to the surface in the absence of 
impurity scattering, will have the same lifetime r = 1/2T. 



Diagonalization of the tight-binding Hamiltonian 

The diagonalization of e is essentially the same as in 
Rcf . |l|. Here we give some details not previously explicitly 
described, and in particular calculate transmission coef- 
ficients across the twin boundary in the superconducting 
state. 

We do not solve the problem for the superconducting 
order parameter near the twin boundary self-consistently, 
but instead substitute a "guess" order parameter into the 
Bogoliubov-de Gennes (BdG) equations 



(16) 



to study the qualitative features of the excitation spec- 
trum. Definitions of nearest-neighbor hopping and pair- 
ing amplitudes are given in Fig. 1. In accordance with the 
orthorhombic symmetry of Cu02 planes in YBCO these 
amplitudes are different along a and h axes: ti^2 = ^(lie) 
and Ai,2 = A(l ±<^). 

In the bulk of each twin, quasiparticles are plane waves 
with excitation energy 



- v/£(k)2 + A(k)2 , 



(17) 



where the normal state quasiparticle dispersion in the 
right (left) twin is given by 



e(k) = —At cos — ;= cos — ^ =F 4ie sm — ^ sm — ^ — ^ 



(18) 



and the mixed symmetry superconducting gap is 



. » ^ . . k'Y-Qi . kiiQi „ . ^ k'Y-Qi 
A(k) - 2Asin^sin^ =F2A(5cos 
V 2 V 2 



ky Oj 



^/2 ^/2 



(19) 

where the first term corresponds to a d-wave component, 
while the second term describes an extended s-wave com- 
ponent (the only s-wave harmonic in our model) which 
admixes with different signs in the two crystal twins. At 
half-filling the maximum amplitudes of the two harmon- 
ics are = A and As = /S.5 < Ad- 

Treating quasiparticle scattering in both the normal 
and the superconducting states one can employ the trans- 
lational invariance along the twin boundary, which leads 



to conservation of the parallel component of the momen- 
tum. For low-energy elastic processes we need to deter- 
mine change in the normal component of the momentum 
of a scattered excitation. From (|l8|), the normal com- 
ponents of the Fermi momenta for incoming and out- 
going electrons with a given parallel component ky are 
— zLk + g, in the right twin, and k'^^ — zLk — q, in 
the left twin (see Fig. 2) . Parameters k and q are defined 
by 



tan 



qa 

71 



e tan 



V2 



ka —fi cos{qa / ^2) 
V2 ~ 4tcos(kya/V2) 



(20) 



The normal components of the Fermi velocities for 
all these states have the same absolute value vfx — 
2\plta smikaj -\/2) cos(fcya/-\/2) / cos(ga/ v^). 

The two-component BdG eigenf unctions on the left 
and right sides, 'tl}f- (x < 0) and ^pf" (x > 0) in the twin 
boundary problem, are still given by a combinations of 
plane waves, therefore among BdG equations ( |l6|) a spe- 
cial care is required only for those with i on the twin 
boundary. Subtracting from the difference equation ( [l^ ) 
for i S TB a uniform part and making Fourier transform 
over y we obtain two conditions: (i) continuity equation 
ipQ — TpQ and (ii) an additional equation "on the deriva- 
tive," which in the leading order in in superconducting 
gap, i.e., neglecting corrections 0{A'^ /ep) to the energy, 
has the following form 

^>e-Wy2 _ ^<,.WV2 ^ Uocos{qa/V2) ^^^^ 

2tcos{kya/V2) 

In the following analysis we will also use the symmetry 
of the twin boundary which is expressed by the relations 
ts-ia-j — tij and Ag-iCTj = —Ay, where a is reflection in 
the twinning plane (not a Pauli matrix) . The BdG equa- 
tions (|l^) are invariant under the combined transforma- 
tion U = TaCT, T3 being the Pauli matrix which acts in the 
particle-hole space. All solutions are, hence, classifled by 
their parity with respect to U . 

In the normal state, quasiparticles move freely across 
the twin boundary for JJq = ^ and acquire a finite scat- 
tering amplitude r — ajii — a) for C/p > 0. 

Let us consider a scattering process below during 
which a quasiparticle approaching the twin boundary 
from the right with the wave-vector is reflected into 
an out-going state with the wave-vector k^ on the same 
side and is transmitted into out-going states with kf 
and k^ on the left side (Fig. 2). We define A_ = A(k>), 
A+ = A(k>), A(k<) = -A+, and A(k<) = -A_, the 
last two relations following from the symmetry of the su- 
perconducting state on the twin boundaryJ3 

There are three cases for the quasiparticle energy E 
to be considered, (i) E < Amin, (ii) '^min < E < A^ax, 
and (iii) Amax < E, where Amin = min(|A_|, |A+|) and 
A,„ax -niax(|A_|,|A+|). 
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1. Bound states, E < Amin 

The bulk states in this energy range are absent and 
the only possibility is states localized in the vicinity of 
the twin boundary. The wave function of a bound state 
in the right twin is given by 



E + iQ 



Q± = JAI-E^, K±=n±/vF. 



(22) 



whereas on the left side is obtained from ipf = 
±T3'i/'?^ for even (odd) eigenstates of the operator U. The 
ratio of the particlc-hole components of the wave- function 
is determined from homogeneous BdG equations, while 
energy E and "reflection" coefficient R have to be found 
from the continuity condition and Eq. (|2l|). Let us con- 
sider these equations in more detail for the odd symmetry 
solutions. From the continuity of the wave function one 
finds 



A_ + RA+ = , 
whereas Eq. (^ij) yields 

E + m^- R{E - in+) = ia[E + in^ 



(23) 

RiE-in+)] , 
(24) 



where a = Uoa/V^vp^. From the real part of the last 
equation we find for the bound state energy 



E = ±- 



2a|A_ 



\/4a^ + [l-R-a^{l + R)]^ 



(25) 



The sign in this equation should be found by substituting 
expression (GS) back into Eq. (p4h. This gives 



\1- R - a^{l + R}\ =T[l-i?-a^(l + ^)] (26) 



and 



R_ _ ^ [a^{l + R) + l~ R] 
Jr] ~ |a2(l + i?) + l-i?| 



(27) 



As follows from Eqs. (|2|) and (|2^) the two cases \R\ < 1 
and \R\ > 1 transform into each other under R ^ 1/R. 
Consequently, we have to consider only — 1 < i? < 1. 

For -1 < i? < or A+A_ > we find from Eq. (H) 
the upper sign if > (1 — i?)/(l + i?) and the lower sign 
in the opposite case. On the other hand Eq. (^ ) gives al- 
ways the lower sign. Therefore, the e nerg ies of the bound 
states of odd symmetry are given by (25) with the minus 
sign as far as the boundary potential is weak and a is 
small. The dispersion of bound state energies disappears 
at Uq = 0. The condition A+A- > is satisfied in the 
vicinity of the gap nodes, e.g., for kyc < ky < kyA and 
for ky > kyD and ky < kyB in Fig. 2. The number of the 
bound states is, thus, proportional to the degree of the 



orthorhombicity of the superconducting gap or As/A^t- 
where As and A^ are amplitudes of s- and d-wave gaps.EJ 
Observation of the zero-energy peak in the local density 
of states near twin boundary would be a direct experi- 
mental evidence of the mixed symmetry gap in YBCO. 

The critical value of the boundary potential Uq, which 
destroys the bound states, is different for states propagat- 
ing normally to the twinning plane and for states prop- 
agating nearly parallel. In the former case a ~ Uq/s, 
where e is geometrical average of the Fermi energy ep 
and the band width. The parameter a is, therefore, quite 
small and cannot exceed the critical value ac = ^- The 
bound states lie near the Fermi level having a weak dis- 
persion given by Eq. ( p5| ) . In the second case of the bound 
states with momenta nearly parallel to the twin bound- 
ary, perpendicular component of the Fermi velocity is 
very small and a is significantly enhanced. Bound states 
near "vertical" nodes in Fig. 2 disappear completely for 
Uo/e > A,/Arf. 

For < i? < 1 or A. A_ < 0, that is, in the "out 
of nodes" region, Eq. ( |27|) yields the upper sign. This 
is compatible with Eq.^2|) if > {1 - R)/{1 + R). 
Arbitrary weak boundary potential can produce bound 
states for the part of the Fermi surface with i? « 1. Total 
number of such bound states is, however, proportional to 
a. As a substitution of R ^ 1 into ( psj) shows, these 
bound states split from the bottom of the continuum 
of bulk quasiparticles, i.e., have energies slightly below 
A(k), which does not vanish in this region of the Fermi 
surface. Therefore, such a possibility being sensitive to 
the value of the boundary potential is less important for 
the thermodynamic of the system. 

Completely analogous consideration of the even sym- 
metry states predicts the upper sign in Eq. ( p5| ) and the 
same condition on the two gaps 



A+A_ > 



(28) 



in order to have bound quasiparticle states on the twin 
boundary. 

Disorder provides an additional mechanism for broad- 
ening the zero-energy peak in the local density of states 
on twin boundary. The critical impurity concentration, 
which smears out completely bound states, is obtained 
by the following arguments. The characteristic space ex- 
tension of bound states is = vf/t^As- In the clean 
limit they can be considered as quasiparticles undergo- 
ing multiple Andreev reflections in the superconducting 
well on the twin boundary. In the presence of impurity 
scattering, the coherence length has to be compared 
to the mean free path I = vft = vf/'^T. For <C ^, a 
quasiparticle is reflected many times before being wiped 
out from the twin boundary region by a scattering on 
an impurity. In this case, though being transformed into 
resonances, the quasibound states are still well defined. 
For > a quasiparticle scatters to the continuum very 
quickly and the twin boundary resonances as well as the 
zero-energy peak disappear. Thus, it is possible to de- 
stroy the zero-energy peak in the local density of states 
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by increasing number of impurities to the critical value 
given by Tc A^. 



2. Extended states with |A_| < E < |A+| 

In this energy range, an excitation, particle or hole, can 
approach the twin boundary from the bulk, but it can- 
not be reflected back into the twin with the wave-vector 
k^, since E < jA^I- The only allowed processes are 
Andreev reflectionMI and Andreev transmissionjd which 
involve particle-hole transformations. The wave function 
of such a quasiparticle is given for a; > by 



-A 



and for a; < 



E. 



B'{ ]e^K '^ 



(29) 



We have neglected in the above formulas small shifts of 
the wave vectors from the Fermi surface for propagating 
states and defined 



E^ 



E. 



E 



Al-E^ . (30) 



Eqs. ( |29| ) describe an incident electron-like quasiparticle 
with momentum k « k^ approaching the twin boundary 
from the right, a reflected hole with wave vector k « 
k^ and transmitted hole with k w k^. There are also 
two damped waves with wave vectors k^ and k^. We 
also took into account the odd symmetry of the gap with 
respect to the reflections in the twinning plane. 

Parameters A, A', B, and B' have to be determined 
from the continuity of the wave function and Eq. (21). 
For zero boundary potential one finds B = B' = so that 
only exponentially decaying wave exists for a; < 0. Hence, 
the transmission coefficient through the twin boundary 
for quasiparticles with |A_| < E < |A+| vanishes. 

For Uq there is a nonzero probability of an incident 
electron at k^ to be transformed into an outgoing hole 
of wave- vector k]^ on the opposite side of the boundary, 
which is a kind of Andreev transmission. Remarkably, 
the boundary barrier Uq > 0, which causes reflection of 
the quasiparticles from the twin boundary in the normal 
state, allows transmission of the low-energy excitations 
in the superconducting state. The particle-hole transmis- 
sion coefficient to the second order in a small parameter 
a is 



Wph 



\B 



'|2 



AXiEl 



,2 \2 



A^A+ + E^E. 



(31) 



where wq — ^o? — 2{U j v f)"^ ■ 
mission coefficient simplifies to 



For < A+ the trans- 



Wpi.^w^iE' - A^_)IE^ 



(32) 



3. Extended states with E > Amax 

The wave function of scattered quasiparticle on the 
left and right sides is given by the same expression ( [29| ) 
except for the decaying factors e*''+^, which are absent 
now, and E+ defined as = E—{E'^ — A^)^/'^ instead of 
Eq. (^0|). The outgoing quasiparticle flow on the left side 
consists now of two parts corresponding to electron and 
hole excitations. Therefore, the transmission coefficient 
has the particle-hole contribution of Eq. ( |3l| ) and the 
particle-particle contribution 



1^ 



A^ 



El _ (Al - El){Al - E^_ 
~W ^ (A_A+ + )2 



(33) 



where the second equality is obtained for a ^ 1. For 
high energies E ^ A,„ax in this case the particle-particle 
transmission coefficient Wpp w 1 and — 0{a'^). 



III. THERMAL RESISTANCE OF TWIN 
BOUNDARY 

In the case of dirty superconductor with strong impu- 
rity scatters, the width F of the low-lying energy levels 
has an important implication on the thermal conductiv- 
ity through the twin boundaries. For F > A^, exci- 
tations with energies approximately equal to F will be 
the most important carriers of heat at low temperatures, 
completely analogous to the case of the bulk thermal con- 
ductivity (see Introduction) . For these excitations w ~ \ 
(if a <C 1) and the flow of the heat will not be much 
affected by the presence of the twin boundaries even in 
the low-temperature limit. 

For F < As, the situation is different. At low temper- 
atures (T ^ F) the heat will again be carried primarily 
by excitations with energies E T and the mechanism 
of the Andreev transmission discussed above produces a 
nonzero flnw of heat across the twin boundary. Following 
Andreev Jl3 the heat current across the boundary is given 
by the expression 



VF = — 2^ E^n^VkxW)^ ■ 



(34) 



Here, the prime attached to the sum over momenta means 
that only excitations with x component of the group ve- 
locity satisfying Vkx < are summed over; the a;-axis is 
the axis normal to the twin boundary. Also, excitations 
with wave vectors near ke and kD in Fig. 2 are neglected 
since their group velocities toward the twin boundary are 
small. The factor 2 accounts for the contribution of two 
spin directions. 

For excitations near kA and kc (Fig. 2), the normal 
projection of the group velocity is 



Vkx 



VFxlEl 



(35) 
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Thus, using Eq. (^2|) for the transmission coefficient we 
have 



(36) 



In integrating over k [on the right hand side of Eq. (34)] 
in the neighborhood of Ica we approximate 



El 



VA^ky (37) 



and \va \ ^ ^o/kp ^ vp. It is easily seen that the effect 
of the factor in brackets in Eq. ( ^6|) , when averaged over k 
as in Eq. (|4|), is to give a factor of the order of unity. For 
example, if the exponent of the factor in square brackets 
were 2 rather than 3, evaluation of the sum over k in 
Eq. ( ^ ) shows that the correct result could be obtained 
by the substitution VkxWji = ^vpwo for Eq. (p^) . Since 
the exponent is 3 not 2, we make a replacement 



VkxWl^ = fvpwo 



(38) 



where / will be close to, but somewhat less than unity. 
Impurity scattering could not change this factor by very 
much. With this substitution and extending the sum over 
k in Eq. (Isl) to all fc, we find 



W 



(39) 



Note, that the quantity in the bracket is £' — fiN, where 
E is the mean energy of the system. 

If the temperature on both sides of the twin boundary 
is the same, the heat current W is balanced by a heat 
current of the same magnitude in the opposite direction. 
Thus, in the presence of a temperature difference across 
the twin boundary, there will be a net heat current given 
by 



ST 



ktb- 



(40) 



where d is an average spacing between twin bound- 
aries. Hence, the thermal conductivity, if limited by twin 
boundary resistance, will have the low-temperature form 



ktb — hCvVpXcS , 



where 



Cv^T 



95 
dT 



V,N 



d_ 

dT 



{E - ^iN) 



and an effective mean free path 
Acff 



-^wod 



(41) 



(42) 



(43) 



has been defined. In writing Eq. (p2|), a weak tempera- 
ture dependence of the chemical potential has been ne- 
glected. 



Note, that since vp and AcS do not depend on T, the 
temperature dependence of the twin boundary thermal 
conductivity will be that of the specific heat. 

In the unitary limit for impurity scattering, the low- 
temperature specific heat at T < F, has the form (|^). 
Thus, the thermal conductivity due to the twin bound- 
aries has the same linear-T dependence at low tempera- 
tures as the bulk thermal conductivity Kbuik- If we now 
consider a sample with one predominant orientation of 
twin boundaries, say parallel to (110), then the thermal 
resistance in the perpendicular direction is additive and 
for the combined thermal conductivity from boundaries 
and impurities we have 



K_L^ = (ktb) ^ + (Kbulk)" 



(44) 



The heat flow in the parallel direction is unaffected by 
twin boundaries and k\\ = Kbuik- Such anisotropy of in- 
plane thermal conductivity disappears for samples with 
equal weights of two types of twins. 

For r ^ r <C As, the specific heat coincides with that 
of a clean d-wave superconductor Cv = 9C(3)A^oT^/Ao. 
Thus both in this case and in the case T > Ag, when 
the scattering can be treated in the Born approximation, 
the twin boundary thermal conductivity varies as T^. 
For the sake of comparison we note that in the unitary 
scattering limit and F ^ T <C Ag the bulk thermal con- 
ductivity varies as T^ (see Introduction), whereas in the 
Born limitliil and for T ^ Aq, Kbuik varies as T. 



IV. CONCLUSIONS 

In the previous work0 on a clean d + s orthorhombic 
superconductor, such as YBCO, we have demonstrated 
the existence of bound zero-energy excitations at the 
twin boundaries and identified an Andreev transmission 
process responsible for heat conduction across the twin 
boundaries at low temperatures. When strong impurity 
scatters are added, a new energy scale F appears, which 
is the relaxation rate of the lowest energy quasiparticles. 
In this paper we show that the twin boundary bound 
states remain well-defined resonances contributing a zero- 
energy peak to the local density of states at the twin 
boundary provided F < A^ , where A^ is an amplitude of 
the s-wave component of the gap. Furthermore, we cal- 
culate the transmission coefficient describing the trans- 
mission of excitations across the twin boundaries, and 
hence evaluate the thermal conductivity of twin bound- 
aries. For F > As, the twin boundaries present resistance 
for the flow of heat due to an extra boundary potential 
only, i.e., the same as in the normal state. For F < A^., on 
the other hand, there are two cases to be considered. If 
T <IC F, the twin boundary thermal conductivity, like the 
bulk thermal conductivity, varies linearly with T, while if 
F ^ T <C As (or if the impurity scatters are weak), the 
twin boundary thermal conductivity varies as T^. The 
best way to identify twin boundary thermal conductivity 
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experimentally would be on a highly twinned sample with 
all twin boundaries parallel to the (110) plane; measure- 
ments of thermal eonductivities in the [110] and the [110] 
directions will then differ due to the thermal resistance of 
the twinning planes. Anisotropy in the basal plane ther- 
mal conductivity, which develops only below T,., would 
be a clear indication of the particle-hole transformation 
in heat flow across twin boundaries. 



FIG. 2. Twin boundary and orthorhombic Fermi surfaces 
in the two twins. The capital letters A, B, ... denote the 

positions of nodes of the d ± s gaps. The sign of A(k) in the 
different regions on the Fermi surfaces is indicated as -|- or — . 
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